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We present for the first time tfie complete matter power spectrum for R" gravity wliicli lias 
been derived from the fourth order scalar perturbation equations. This leads to the discovery of a 
characteristic signature of fourth order gravity in the matter power spectrum, the details of which 
have not seen before in other studies in this area and therefore provides a crucial test for fourth 
order gravity on cosmological scales. 



Ever since the Concordance model yy] was proposed as 
the best fit to all available cosmological data sets, there 
have been many attempts to understand the nature of 
Dark Energy. However, despite enormous effort over the 
past few years, this problem remains one of the greatest 
puzzles in contemporary physics. One of the theoretical 
proposals that has received a considerable amount of at- 
tention recently, is that Dark Energy has a geometrical 
origin. This idea has been driven by the fact that modi- 
fications to General Relativity appear in the low energy 
limit of many fundamental schemes and that these 

modifications lead naturally to cosmologies which admit 
a Dark Energy like era [4j without the introduction of 
any additional cosmological fields. Most of the work on 
this idea has focused on fourth order gravity, in which the 
standard Hilbert-Einstein action is modified with terms 
that are at most of order four in the metric tensor. The 
features of fourth order gravity have been analyzed with 
different techniques and all these studies suggest that 
these cosmologies can give rise to a phase of accelerated 
expansion, which is considered to be an important foot- 
print of Dark Energy. 

The work described above has largely focused on the 
dynamics of homogeneous cosmologies which have the 
standard Robertson- Walker geometry and are therefore 
also isotropic. Although these results have many of the 
desirable features that we are looking for, such as a mat- 
ter dominated epoch and late-time acceleration, there 
are still some key issues that need to be addressed be- 
fore one could claim to have a cosmological description 
which is able to compete with the standard ACDM cos- 
mology. The calculation and analysis of the evolution 
of linear perturbations and their comparison with ob- 
servations is clearly among the most important of these 
open problems. Over the past year this problem has been 
studied by several authors using the metric approach to 
perturbations, by either considering different ways of pa- 
rameterizing the non-Einstein modifications of gravity or 
by simplifying the underlying fourth-order perturbation 
equations using a quasi-static approximation [6,]. 



In what follows, we demonstrate that considerable 
progress can be made to this problem by using the 1+3 
covariant approach to cosmological perturbations 7] . Us- 
ing a specific recasting of the field equations (based on the 
Ricci and Bianchi identities), the development of equa- 
tions describing cosmological perturbations in theories of 
gravity characterized by an action which is a general an- 
alytic function of the Ricci scalar f{R), becomes both 
transparent and straightforward, allowing for the exact 
integration of the perturbation equations without mak- 
ing any approximations. In order to easily discuss the key 
features of the perturbation dynamics and the associated 
power spectrum, we focus on i?"-gravity. This theory is 
characterized by the action L = \J—g [xR" + ^m] and is 
the simplest possible example of fourth order gravity. 

Before we can discuss the evolution of density pertur- 
bations, a suitable background cosmology must first be 
found. In d, Q , the complete dynamics of homogeneous 
and isotropic cosmologies were studied in detail using 
the dynamical system approach (see and references 
therein). It was found that in i?" gravity, it is possible 
to have a transient matter-dominated decelerated expan- 
sion phase, followed by a smooth transition to a Dark 
Energy like era which drives the cosmological accelera- 
tion. The first phase, characterized by the baratropic 
equation of state parameter w, has an expansion history 
determined by a scale-factor a = ^211/3(1-1-10) .^^jj^gj-g .^^g j-g_ 
strict ourselves to n > for this background as negative 
values of n would represent a contracting model. This 
solution provides exactly the setting during which struc- 
ture formation can take place and is therefore an ideal 
background solution for our study of density perturba- 
tions. 

Scalar perturbations, which describe density perturba- 
tions may be extracted from any first order tensor Tab 
orthogonal to by using a local decomposition [ll| , so 
that repeated application of the operator Va = h^^b on 
Tab extracts the scalar part of the perturbation variables. 
In this way we can define the following scalar quantities 
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where A™, Z respectively represent the fluctuations in 
the matter energy density and expansion O, and TZ, 
3? determine the fluctuations in the Ricci scalar R and 
its momentum R. This set of variables completely char- 
acterizes the evolution of density perturbations. Then, 
using eigenfunctions of the spatial Laplace-Beltrami op- 
erator defined in V^Q — —-^Q, where k — 2nS/X is 
the wave number and Q = 0, we can expand every first 
order quantity in the above equations: 



X(t,x)-^xW(t) QW(x), 



(2) 



where ^ stands for both a summation over a discrete 
index or an integration over a continuous one. In this 
way, it is straightforward, although lengthy, to derive a 
pair of second order equations describing the fc*'* mode 
for density perturbations in /(i?) gravity. They are: 
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where /' = df{R)/dR, the quantities fiR, pR are the en- 
ergy density and pressure of the curvature fluid defined 
in [ij] and if = 0, —1 is the usual spatial curvature 
scalar . It is easy to see that for the f{R) = R case, these 
equations reduce to the standard equations for the evo- 
lution of the scalar perturbations in General Relativity. 

Already on super-Hubble scales, k/aH <C 1, a number 
of important features are found which allows one to dif- 
ferentiate ^ from their General Relativity counterparts 
p^ . Firstly, it is clear that the evolution of density per- 
turbations is determined by a fourth order differential 
equation rather than a second order one. This implies 
that the evolution of the density fluctuations contains, in 
general, four modes rather that two and can give rise to 
a more complex evolution than the one of General Rela- 
tivity (GR). Secondly, the perturbations are found to de- 
pend on the scale for any equation of state for standard 
matter (while in General Relativity the evolution of the 
dust perturbations are scale-invariant). This means that 
even for dust, the evolution of super-horizon and sub- 
horizon perturbations are different. Thirdly, it is found 



that the growth of large density fluctuations can occur 
also in backgrounds in which the expansion rate is in- 
creasing in time (see flgure 1). This is in striking contrast 
with what one flnds in General Relativity and would lead 
to a time-varying gravitational potential, putting tight 
constraints on the Integrated Sachs- Wolfe effect for these 
models. 

Let us now turn to the case of a general wave mode 
k. One of the most instructive way of understanding 
the details of the evolution of density perturbations for 
a general k is to compute the matter power spectrum 
P(fc), deflned by the relation [l3| (A„(ki)A„(k2)) = 
P(fci)(5(ki -l-k2), where k^ are two wavevectors character- 
izing two Fourier components of the solutions of ([3]) and 
P(ki) — P{ki) because of isotropy in the distribution of 
the perturbations. This quantity tells us how the fluctu- 
ations of matter depend on the wavenumber at a specific 
time and carries information about the amplitude of the 
perturbations (but not on their spatial structure). In 
General Relativity, the power spectrum on large scales is 
constant, while on small scales it is suppressed in compar- 
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FIG. 1: Plot against n of the real part of the long wavelength 
modes for i?"-gravity in the dust case (blue, red green and 
yellow lines) together with the GR modes (red and purple 
line). Note that there is at least one growing mode for any 
value of n. This means that even in cases where the expan- 
sion rate is accelerating, i.e., n > jiT^' large-scale density 
perturbations grow. 



ison with the large scales (i.e., modes which entered the 
horizon during the radiation era) [14] . In the case of pure 
dust in General Relativity the power spectrum is scale in- 
variant. Substituting the details of the background, the 
values of the parameter n, the barotropic factor w, the 
spatial curvature index K and the wavenumber k into ([3]) 
one is able to obtain P{k) numerically. 

The k-structure of equations ^ suggest that in fourth 
order gravity there exist at least three different growth 
regimes of the perturbations. This is confirmed by our 
results (see figure 2). In particular, in the case of dust 
we have three regimes for any values of the remaining 
parameters: (i) on very large scales the spectrum it is 
like what one finds for General Relativity, i.e., scale in- 
variant; (ii) as k becomes bigger the scale invariance is 
broken and oscillations in the spectrum appear; (iii) for 
even larger k the spectrum becomes again scale invariant. 
However, on these scales the spectrum can contain either 
an excess or deficit of power depending on the value of 
n. In particular for n « 1"*" small scales have more power 
than large scales, but, as one moves towards larger values 
of n, the small scale modes are suppressed. 

The features of the spectrum that we have presented 
above can be best interpreted by comparing the sys- 
tem ([3]) which produced it with the equations for the 
evolution of scalar perturbations for two interacting flu- 
ids in General Relativity [lB|. Immediately one notices 
that they have the same structure, i.e., there are friction 
and source terms due to the interaction and the gravi- 
tation of the two effective fluids. It is then natural to 
ask ourselves if this analogy can be useful to better un- 
derstand the physics of these models. The answer is af- 
flrmative. First of all a more correct way to draw this 
analogy would be to write the system of equations for 
Am and fluctuations in the energy density of the curva- 
ture fluid An = S'^'V^ fin/ and analyze their structure 
rather than using the ones above. On very large and on 




FIG. 2: Plot of the power spectrum at r = 1 for i?"-gravity 
and n > 1. Note that the spectrum is composed of three 
parts corresponding to three different evolution regimes for 
the perturbations. 
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FIG. 3: Evolution the Power spectrum for /^"-gravity for 
n = 1.4. The spectrum has been normalized in such a way 
that the curves coincide at large scales. Note how, as time 
passes, small scale perturbations are dissipated and oscillation 
appear. 



very small scales, the coefficients of the (A^, Aj^) sys- 
tem become independent of k, so that the evolution of 
the perturbations does not change as a function of scale 
and the power spectrum is consequently scale invariant. 
On intermediate scales the interaction between the two 
fiuids is maximized and the curvature fluid acts as a rel- 
ativistic component whose pressure is responsible for the 
oscillations and the dissipation of the small scale pertur- 
bations in the same way in which the photons operate in 
a baryon-photon system. This suggests the following in- 
teresting interpretation for the perturbation variables TZ 
and 3?. These quantities can be interpreted as represent- 
ing the modes associated with the contribution of the ad- 
ditional scalar degree of freedom typical of /(i?)-gravity. 
In this sense the spectrum can be explained physically 
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as a consequence of the interaction between these scalar 
modes and standard matter. The result is a consider- 
able loss of power for a relatively small variation of the 
parameter n. For example, in the case n — 1.4 the differ- 
ence in power between the two scale invariant parts of the 
spectrum for n = 1.1 is of one order of magnitude while 
for n = 1.6 is about ten orders of magnitude. It should 
be noted that existing analysis [l^ of this model require 
n = 3.5 in order for predictions to be consistent with 
measurements of rotation curves of low surface bright- 
ness galaxies and SNe la. Given the huge drop in power 
at small scales in the power spectrum for n — S.b one 
expects that this model could be easily ruled out. 

Further information on the dynamics of the matter per- 
turbations can be obtained examining the time evolution 
of the power spectrum. In figure [3] we give the power 
spectrum for n = 1.4 at different times. One can see 
that, as the universe expands, the small scale part of the 
spectrum is more and more suppressed and oscillations 
start to form, suggesting that in this model small scale 
perturbations tend to be dissipated in time. On the other 
hand on large scales they do not evolve, which might ap- 
pear in contrast with what mentioned above. However 
this is a byproduct of the normalization: for clarity we 



have normalized the spectrum in such a way that every 
curve has the same power in long wavelength limit. A 
more in depth discussion of the features presented above 
can be found in [17]. 

Probably the most important consequence of the form 
of the spectrum presented above is the fact that the ef- 
fect of these type of fourth order corrections is evident 
only for a special range of scales, while the rest of the 
spectrum has the same k dependence of GR (but differ- 
ent amplitude). This implies that we have a spectrum 
that both satisfies the requirement for scale invariance 
and has distinct features that one could in principle de- 
tect, by combining future Cosmic Microwave Background 
(CMB) and large scale surveys (LSS) [li,[i|. 
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